Perturbation-induced radiation by the Ablowitz-Ladik soliton 
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An efficient formalism is elaborated to analytically describe dynamics of the Ablowitz-Ladik soli- 
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I. INTRODUCTION 

Dynamics of discrete solitons (intrinsic localized 
modes) in nonlinear lattices has become a topic of in- 
tense research summarized in a number of excellent re- 
views 0. Propagation properties of waves arising as a 
result of the interplay of nonlinearity with lattice dis- 
creteness can be quite distinct from those inherent in 
continuous nonlinear systems and hold much promise for 
applications in various physical, biological and techno- 
logical problems. Examples are energy localization and 
transfer in systems of nonlinear oscillators 0], propaga- 
tion of self-trapped beams in arrays of coupled nonlinear 
optical waveguides nonlinear charge and excitation 

transport in biological macromolecules [M^i local denat- 
uration of DNA double helix 0, dynamics of localized 
excitations in arrays of coupled Josephson junctions 
propagation of optical spatial solitons in nonlinear pho- 
tonic crystals |2l and in diffraction-managed waveguide 
systems [Tol . creating discrete solitons in Bose-Einstein 
condensate Recently it was proposed 1^ to use 

discrete solitons in two-dimensional networks of nonlin- 
ear waveguides to realize functional operations such as 
blocking, routing, logic functions, and time gating. 

Most of the above phenomena are modelled by the 
discrete nonlinear Schrodinger (DNLS) equation or, in a 
more general setting, by the discrete self-trapping equa- 
tion I2I . Recent developments in the study of the DNLS 
equation are reviewed in Refs [T^, 14]. However, the stan- 
dard DNLS equation is nonintegrable jl^ 0| and does 
not exhibit exact soliton solutions, though it can be de- 
rived as a discretization of the integrable continuous NLS 
equation. Hence, numerical methods are generally used 
to investigate nonlinear lattice dynamics on the basis of 
the DNLS equation. 

On the other hand, there exists the integrable dis- 
cretization of the NLS equation - the Ablowitz-Ladik 
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(AL) equation fl^ which has exact soliton solutions and 
admits the complete description in the framework of the 
inverse spectral method. Moreover, Konotop et al. 
and Cai et al. 19] proved integrability of the inhomoge- 
ncous AL system in an external electric field of a partic- 
ular form. Being unique from the mathematical point of 
view, the AL equation is less applicable in physics than 
the DNLS equation. Salerno 20] introduced an equation 
that interpolates between the DNLS and AL equations 
and permits studying (as a rule, numerically) the role 
of integrable and non-integrable contributions to lattice 
properties . The AL-DNLS system with an impurity 
was investigated by Hennig et al. ^| 

A different point of view on the interrelation between 
the AL and DNLS equations was posed in Refs. [23,13, 
I25I ]. In a definite region of parameters the DNLS equa- 
tion can be treated as a perturbed version of the AL 
equation. When a perturbation is small, the discrete 
soliton perturbation theory can be successfully applied 
to analytically describe localized excitations in a system 
governed by the DNLS equation. Such an approach was 
developed in Refs. 0, US in the framework of the 
adiabatic approximation, when a perturbation-induced 
radiation is ignored and a perturbation manifests itself as 
a slow evolution of initially constant AL soliton param- 
eters. The evolution equations for the parameters were 
derived by Vakhnenko and Gaididei 23]. Stability as- 
pects of Hamiltonian perturbations for the AL equation 
were discussed by Kapitula and Kevrekidis 26] . Recently 
the perturbative method to study the AL soliton dynam- 
ics was used in Ref. 27] in relation to energy transport 
in a-helical proteins and in Ref. [2^ for the soliton in a 
random medium. Besides, AbduUaev et al. |29l ] proved 
the existence of discrete autosolitons in the AL model 
with linear and quintic damping, cubic amplification and 
complex filtering (the discrete complex Ginzburg-Landau 
model). Exact solutions of this model for certain rela- 
tions between parameters are given in Ref. [3fl|. 

It is well known that a perturbation of the soliton is 
also accompanied by radiation of small-amplitude disper- 
sive linear waves (or shape distortion) |3l| . and a com- 
plete description of the perturbed soliton dynamics ne- 
cessitates accounting for both the soliton parameter evo- 



2 



lution and the radiation effects. Therefore, the main goal 
of the present paper is to develop a corresponding (rela- 
tively simple) formalism and to extend, as far as possible, 
the applicability of analytical methods in studying near- 
integrable nonlinear discrete systems. It should be noted 
in this connection that Konotop et al. 32] derived by 
means of the Gel'fand-Levitan-like summation equations 
the evolution equation for the reflection coefficient in the 
case of the inhomogeneous AL model but without using 
it for specific calculations. An estimation of radiative 
corrections to the AL soliton subjected to the stochas- 
tic perturbation was outlined in the important paper by 
Garnier on the basis of conserved quantities. 

Our approach utilizes the Riemann-Hilbert (RH) prob- 
lem "33^ . The application of the RH problem to perturbed 
nonlinear equations was initiated by Kivshar on an 
example of the Landau-Lifshitz equation. A purely al- 
gebraical calculation of higher-order corrections to the 
perturbed NLS soliton and of radiation effects for a soli- 
ton in a doped fiber was performed on the basis of the RH 
problem in Ref. . Such an approach has been proved 
to be efficient for a wide class of continuous perturbed 
nonlinear equations, including multicomponent ones |36j ]. 

This paper gives a self-contained exposition of the AL 
soliton perturbation theory. In Sec. ^we fix preliminary 
facts concerning the AL spectral problem which are used 
in Sec. IIIII to formulate the RH problem. In Sec. IIVI we 
describe a procedure to solve the RH problem with zeros 
and obtain immediately the AL soliton solution in Sec. 
IVl We stress that calculations within the RH problem 
do not use discrete analogs of the Gel'fand-Levitan inte- 
gral equations. Section|^is devoted to derivation of the 
evolution equations for the RH problem data associated 
with the AL soliton parameters. These equations exactly 
account for the perturbation and serve in the subsequent 
Sections as the generating equations for the perturbative 
expansion. Sec. IVIII contains brief exposition of the adi- 
abatic approximation, whereas Sec. IVIIII represents the 
main result of the paper - derivation of formulas for cal- 
culating radiative corrections from the continuous part 
of the RH problem data. In Sec. IIXI we illustrate the 
formalism by the examples of linear damping and quintic 
perturbations. Appendices contain some technical details 
of the applications of the RH problem. 

II. THE ABLOWITZ-LADIK SPECTRAL 
PROBLEM 

A. Jost solutions and asymptotics 

Integrable discretized NLS equation (AL equation) 
iu„f + u„+i+u„_i-2u„ + |u„P(u„+i+u„-i) = (2.1) 

for a scalar complex function u which depends on discrete 
variable n, — cx) < n < 00, and time t admits the Lax 
representation with the AL spectral problem 17] 

J{n + l)^{E + Qn)J{n)E-^, (2.2) 



/Om„\ ( z Q \ 

^" - 1^ -< j ' ^ " 1^ z-i j ' 

and the evolutionary equation (subscript t means time 
derivative) 

J^{n) = V{n)J{n) - J{n)VL{z), (2.3) 
Vin) = i( ~ ^"'""+0 + r!, 

m = ^{z-z-'^as. 

It means that Eq. (|2.1|l arises as a compatibility condition 
for Eqs. H2.2|l and H2.3|l . Here z is a constant spectral pa- 
rameter and the star stands for the complex conjugation. 
The spectral problem in the form (|2.2|l differs slightly 
from the usual one |l7j and permits introducing matrix 
Jost solutions J±{n) of Eq. (|2.2|l with the unit asymp- 
totics, J±(n) 1 as n — ^ ±00. J±{n) solve Eq. H2.3|l as 
well. The scattering matrix >5'(2:) defined by 

J_(n) = J+{n)E''S{z)E-"- (2.4) 

has the structure 

The Jost solutions obey the conjugation condition 

Jl{n,z)^v±{n)j^\n,z), (2.5) 

where z — 1/z*, 'f means the Hermitean conjugation 
and 

00 n— 1 

«+N = ^^^^ w_(n)= ]j pi, pi = i+\ui\^. 

l—n I — — 

We also obtain that detJ±(n, z) — v±{n), detS' = -y, 
where v — Il^-oo ^^^^ evidently v+{n)v — v^{n). 
From Eqs. 12.4|l and (|2.5(l we obtain S^z) = vS~'^{z) 
which gives a+(z) — a_(z), 6+(z) — b^{z). 

The AL spectral problem 12.2|l obeys the important 
symmetry ('T-'-parity'): if J{n,z) is a solution, then 

VJ{n, z) = a3J{n, ~z)(T3 (2.6) 

is a solution, too. It follows from Eq. 12.6|l that diagonal 
elements of J(n, z) are even functions of z, while off- 
diagonal entries are odd functions. The same symmetry 
is valid for the Jost solutions and the matrix 5, the latter 
means a±{z) ~ a±(— ^), b±{z) = —b±{—z). 

Now we consider asymptotic behavior of the solution 
J(n, z) for z — > 00. Let 

J{n,z) = j(°)(n) + z-ij(i)(n) + 0(z-2), z ^ 00. 

Inserting this expansion into Eq. (|2.2() gives 

J(")(n + l)=(j J(«)(n), (2.7) 
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while the potential w„ is retrieved as 

- _ 7(1) / 7-(0) 
"^12 / "^22 ■ 



det*~ 



v_^{n)a-{z), 



T+S = T_ 



(2.8) 



Note that the asymptotics (|2.7|l is consistent with the 
■P-parity property (|2.6|l . Similar results hold for z ^ 0, 
when J{n, z) — J(-o) (n) + zjfi^ (n) + 0{z'^): 

J(o){n + l) = 0" 1 ^ J(o)in). 



B. Analyticity 

Let C± are domains in the complex z-plane lying out- 
side (+) and inside (— ) the unit circle \z\ = 1. It fol- 
lows from the spectral problem (|2.2|l that the first col- 
umn {n, z) of the Jost function J_ and the second 

rrpi 

one J\_ {n, z) of J+ are analytical in C+ (and continuous 
for z —>■ 1). Hence, the matrix function 

M/+(n,z)= (jW,ji")(n,^) 

is a solution of the spectral problem H2.2|) and analyti- 
cal as a whole in C+. We obtain from the conjugation 
formula Eq. I|2.5|l that the rows {J-)'^^ and iJ+)^2] ^^"^ 
analytical in C_. As a result, the matrix function 



is analytical as a whole in C_ and solves the adjoint spec- 
tral problem. 

Analytical solutions can be expressed in terms of the 
Jost functions. Indeed, 



as well as 



a+ 



J+E-S+E-"^ S+{z) - 5+ 1 , (2.9) 



1 b_/v 
a+/v 



= J^E^'S-E-'^, S- : 
= SS — 

It follows from these formulas that 

det*+(n, z) = w+(n)a+(z). (2.10) 
In the same way we obtain 

^-1 ^ E"T+E-"J^'^ = E^T-E-'^JZ^, (2.11) 



r+ = 



1 
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Asymptotic behavior of analytical solutions is derived 
directly from that of the Jost functions and Eqs. H2.9|l 
and lEmj: 



*+(n,z) 
^Z^{n, z) 



1 

v+{n) 

vZ^{n) 







1 



cx), (2.12) 



0. 



Hence, det 5*+ v^{n) as z — > 00 which gives from 
Eq. I|2.1()(l a+{z) ^ 1 as z — > 00. Similarly, a_(z) 1 
as z — > 0. The conjugation formula for the analytical 
solutions follows from Eq. I|2.5|l : 

V-{n) 
v+{n) 

(2.13) 



^fUn, z) = B(n)*li(n, z), B[n) 



III. MATRIX RIEMANN-HILBERT PROBLEM 

Having matrix functions '^j^ and '^Z^ which are ana- 
lytical in two complementary domains C+ and C_ of the 
z-plane and continuous on the contour |z| = 1, we can 
pose the matrix Riemann-Hilbert (RH) problem 

*:^(n,z)*+(n,z) =i;"G(z)£;-", |z| = 1 (3.1) 

as a problem of analytical factorization of the matrix 
function G(z) defined on the unit circle |z| = 1. It follows 
from Eqs. 1)2. 9|l and H2.11|l that the matrix G has the form 



1 b^/v 
b+ 1 



(3.2) 



The normalization of the RH problem is given by 
Eq. l(TT^ . 

The RH problem (|3.1() has a non-canonical normaliza- 
tion depending on the potential w„. It has been proved 
[3^ that it is possible to reformulate the AL spectral 
problem H2.2|l so as to arrive at the RH problem with the 
canonical normalization and to give a Hamiltonian for- 
mulation with the canonical Poisson brackets. However, 
the above canonicity is achieved at the cost of nonlinear 
dependence of the spectral problem on the potential. Be- 
ing useful for treating non-perturbative AL equation and 
its integrable generalizations, such an approach seems to 
be of less value for the case of perturbations. 

In general, the matrices 4*+ and have zeros in 
some points Zj and Zk in their regions of analyticity. 



i.e., det\E'+(zj) — 0, zj G C+, j = l,2,...,Af+, and 
det*l\zfc) = 0, Zfe e C_, fc = 1,2,...,7V_. We 
suppose that all zeros are simple and of finite number 
with iV+ = 7V_ = N (in other words, we have zero- 
index RH problem). In virtue of the P-parity, zeros 
appear in pairs as ±Zj and ±Zfc. Taking into account 
Eqs. (|2.1U|) and H2.11|l . we conclude that zeros of ^f^. and 
^I'Z^ are given by zeros of the scattering matrix elements: 



a+(±Zj) = and a_(±Zfe) 



0. 
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IV. REGULARIZATION OF THE 
RIEMANN-HILBERT PROBLEM 

We will solve the RH problem H3.1|l with zeros by 
means of its regularization. This procedure consists in 
extracting from rational factors which are responsi- 
ble for the existence of zeros. Indeed, near the point Zj 
we have det ^{z) ~ {z — Zj). Let us introduce a rational 
matrix function SJ^ = + — Zj)(z — Zj)^^Pj, where Pj 
is a rank 1 projector, P'j = Pj. Because Pj — diag(l,0) 
in an appropriate basis, we obtain det = (z — Zj)[z — 
Zj)~^. Therefore, the product 5'-|.(z)SJ^(z) is regular 
in the point Zj (its determinant is nonzero in this point). 
Regularization of the zero —Zj is given by a rational func- 
tion SI^ = 1 — {zj — Zj){z + Zj)~^P^j. As a result, the 
matrix function z) — 5'_(.(n, is regular 

in the points . In the same manner we regularize the 
matrix '^Z^{n^ z) with zeros in ztz^. Namely, the matrix 



is regular in the points ±Zfc and 
Zk - Zk 



^k 



a 



Z - Zk 



„ _ Zk — Zk „ 

Pk, E^k = ^ + ^—^P-k- 

Z + Zk 



Regularizing all iN zeros of the RH problem ()3.1|) . we 
represent the functions as a product 



of the rational matrix function 



(4.1) 



r(7i, z) = ^-n^nE-{n-i)^n-i ■ ■ ■ ^-iSi (4.2) 

and the holomorphic matrix functions ip±{n,z) which 
solve the regular RH problem (i.e., without zeros): 

ipZ\n, z)iP+{n, z) = r(n,z)£;"G(z)i;-"r-i(n,z). 

The appearance of a simple zero Zj of the matrix '5+ 
means that there exists an eigenvector \ j) with zero eigen- 
value. 



«'+(n,z,)|j) =0. 



(4.4) 



Taking the Hermitean conjugation of this equality with 
account of the conjugation property (|2.13|) . we obtain 
(j|i35'I^(n, Zj) = with (j\ — Therefore, the pro- 
jector Pj can be naturally defined as 

\j){j\B 



P. 



{j\B\j) 



(4.5) 



For the zero —Zj we have ^'+(?t., ^Zj)\ — j) = 0. In virtue 
of the 'P-parity, both vectors \ j) and | — j) are interrelated, 
\-j) = <^3\j), and therefore P-j = (J3Pj(T3- 



For practical purposes, it is more convenient to decom- 
pose the products H4.2II into simple fractions. Following 
Refs. USE^, we obtain 

2N 

r(n,z) - 1 - ^ —^\y,){D-^),k{yk\B, (4.6) 

Z Zk 



T-^ln-z) 



2N ^ 

1+ E ^—\yj){D-^)Ayk\B 



z — z^ 

j,k=l ^ 



with new vectors \yj), where zeros are enumerated as 
zi, — zi, Z2, — Z2, . . . , zat, — zat, (and similarly for ±Zfc), 
whereas matrix elements Dkj are given by 



Dkj = (yk\ — ^^1%)- 

Zj Zk 



(4.7) 



It is seen from Eq. (|4.t)|l that the asymptotic expansion 
for r(n, z) has the form 

r(n,z) = 1 + z-^T'^^\n) + 0{z-^). 

Because = + z^i*'^^^ + 0{z-^) = ^+(1 + 

2;~ir(^^ + 0{z^^), we can choose a z-independent func- 
tion V-'+ as a solution of the regular RH problem (|4.3|l : 





-in) 



(4.8) 



where the last equality follows from Eq. (|2.12|) . There- 
fore, in accordance with Eqs. H2.8I) and (|4.1I) . the solution 
Unit) of the AL equation can be retrieved from the solu- 
tion of the RH problem as 



Unit) = — lim 



iz^+) 



12 



(1) 



r 



(1) 

12 



+22 



(0) 
+22 



v+{n) ' 



(4.9) 



The matrix T is mainly determined by the vector \yj). 
Now we derive a coordinate dependence of the vector. It 
follows from the spectral problem that 

'^+{n + l,Zj)\y^,n+l,t} = 

= {E{zj) + Qn) *+(n, Zj)E-^{zj)\yj,n+ l,t}. 

Hence, we can pose E^^{zj)\yj,n + l,i) = |2/j,n,f}, or 

|y„n,i)=£;"(zj)|y„i), (4.10) 

where the vector \yj,t) does not depend on n. Simi- 
larly, it follows from Eqs. H2.3|l and (|4.4|) that n, t)t = 
57(zj)|j/j, n, i). Therefore, the coordinate dependence of 
\yj,n, t) is given as 

\yj,n,t) = E"{zj)e^^'''^*\p), \p) = const. (4.11) 

Finally, we find from the identity det '^+{zj,n,t) — 
that zeros Zj do not depend on n and t. Zeros ±Zj, ±Zj 
and vectors \yj,n,t) comprise the discrete part of the 
RH problem data, while the functions b±{z) entering the 
matrix G (|3.2() are responsible for the continuous data 
with the dependence on t of the form 



Gt = [n,G]. 



(4.12) 
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FIG. 1: Typical arrangement of zeros corresponding to a sin- 
gle soliton. 



V. SOLITON SOLUTION 

In what follows we will not dwell on evident formulas 
for the general case of 4Af zeros. Instead we will give a 
detailed account of the case of four zeros corresponding to 
a soliton (Fig. 1). Hence, after the regularization of the 
matrix RH problem with zeros ±zi and ±zi, we arrive 
at the regular RH problem (|4.3|l . Solitons are associated 
with the discrete part of the RH data, while the continu- 
ous data are now trivial (G = 1). Hence, the solutions of 
the regular RH problem are written in accordance with 
Eq. 1^ as 

^+(n)=^_(n)=(j J^)). (5.1) 

It is possible to express v^(n) (and w_(n)) through 
r(n, z ~ 0). Indeed, because now 5*+ = we find 
from Eq. (|2.12|) '^^ — > diag(w_ (n), 1) as z ^ which 
gives from = t/i+E and Eqs. (|4.11|) and %l.Vi\ : 

E(n,0) = diag(i;_(n),^;+i(n)), (5.2) 
B = diag(En(n,0),E22(n,0)-i). 

Thus, the reconstruction formula H4.9II for solitons is writ- 
ten more conveniently as 



u„(t) = -Ef2^(n)E22(n,0). 



(5.3) 



Hence, it is the matrix E(n, z) that determines the soli- 
ton solution. For simplicity, we will hereafter denote the 
vector |yi,n,t) as \n) with |j/_i,n, i) = CTsIn). Denoting 
Zl = exp[(l/2)(/j,-t-iA:)] and (j>i/p2) = exp{a + i(p), where 
Pi and p2 are components of the constant vector \p), we 
find from (I4.11|) the vector |n) explicitly: 



In) 



62 



g5(2;„+iV,.) 
g-5(2;„+iV>.) 



(5.4) 



Here x„ = /in — 2tsinh/isinfc + a, — kn 

2i(cosh/iCos/c — 1) + Lp. 



As regards the matrix E, it follows from Eq. (|4.6|l with 

= 1, Z2 = — Zl and Z2 = — zi that 



E(n,z) = 1 
1 



(5.5) 



Z — Zl 

1 



[\n){D-%^){7i\B + aM{D-%MB\ 
[\n){D-%2{n\B<yi + a:i\n){D-%2{n\B<J:i\ 



Z + Zl 

Calculating then matrix elements Di^j (|4.7|l with 
detE(n,0) = exp(2/i), we obtain from Eq. (|5.5|l 

/ N sinh « ~ , , sinh u ~ , ^ 

E(n,z)^l- ^ F„(n)- ^. ,:, f+(n), (5.6) 



2(z-zi) 



2(z + zi) 



^ 1 , s sinh u. ^ , , sinh u „ , , 



2(z + zi) 



where 



cxp|^/j.(n a:) + -|A;j cxp[i/c(n — a;)+'iQ — /x] 

p I cosh /x(n— 1— a:) cosh /x(n— 1— x) 

cosh fi{n—x) cosh fi{n—x) 

cx.p^fi{n— ^—x)-\-^k^ cKp[ik{n—x)-\-ia. — fi\ 

/ \ I cosh fi{n—x) cosh fi{n—l—x) 

' ^ cxp[- ik(n-l-x)-ia+^j] exp[-^(n- 1 -a:)+ f fc] 

cosh fi{n—x) cosh fi{n—l—x) 

F+{n) = -a3F_in)a3, F+{n) = ~a:iF_{n)a^. (5.7) 
Here 

, , sinh/x 3 , , 

xyt) — 2t smk + xo, xq ^ -, (5.8) 

ji /i 2 

k 

a{t) — 2t(cosh/icosA: H — sinh /i sin fc — 1) + ag, 
ak 

Uq = if k. 

M 

As a result, we obtain from Eq. (|5.3fl the AL soliton so- 
lution fl7| : 



'Un(i) = exp[ifc(n — x) -\- ia] 



sinh /i 



(5.9) 



cosh/i(n — a;) 

Here and in what follows we write for simplicity 
cosh[/i(n — x)] as cosh/i(n — x). The AL soliton de- 
pends on four constant parameters /i, fc, Xq and ag 
which determine soliton mass 2/i, its group velocity 
Vgr = 2(sinh/x//i) sinfc, soliton maximum position x{t) 
and phase a{t). 

It should be noted for later use that in the presence 
of a perturbation Eqs. H5.8I) are modified due to possible 
perturbation-induced evolution of the soliton parameters: 

2 /■* 

x(t) = - / s\n\iiji{t')s\nk{t')dt' + xo{t), (5.10) 
M J 

a{t) = J [cosh p{t') cos k{t') - l]dt' 
k /■* 

-I- 2- / sinh fi{t') smk{t')dt' + ao{t). 
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VI. PERTURBATION-INDUCED EVOLUTION 
OF THE RH DATA: EXACT RESULTS 

Having formulated the basic ingredients of the RH ap- 
proach to the AL system, we now proceed to the consid- 
eration of the perturbed AL equation 

iu„t + Un+i + u„_i - 2u„ -I- |u„P(m„+i -I- u„-i) = er„. 

(6.1) 

The small parameter e characterizes the perturbation 
amplitude and r„ describes the functional form of the 
perturbation. To find corrections to the soliton caused by 
a perturbation, we first derive the corresponding evolu- 
tion of the RH data. In order to distinguish between the 
'integrable' and 'perturbative' contributions to the evo- 
lution equations, we will assign the variational derivative 
d/5t to the latter. For example, we write iSun/5t — er„, 
as follows from Eq. (|6.1I) . or, in matrix form. 



Rn — 



r„ 



(6.2) 



A. Continuous data 

Consider the spectral problem (|2.2ll . The perturbation 
causes a variation SQn of the potential which in turn 
leads to a variation of the Jost solutions. It follows from 
Eq. H2.2|l that these variations are written in the form 

E~"Jl\n)6J-in)E'' (6.3) 
= E-^'+^\jZ\l + l)SQiJ-{l)E\ 

I — — CO 

E-'\J^\n)SJ+in)E'' 

OO 

= -J2e-^'^'^J+\1 + 1)6QiJ+{1)E\ 

l—n 

where SQi = {6Qi/5t)6t and we have used dJ±{n) ^ as 
n ±oo. Hence, due to the definition H2.4|l . we obtain 
from Eq. H6.3|) a variation of the scattering matrix: 

— = -ieS+T+{z)SZ^ = -ieT^^T^{z)T^. (6.4) 

Here the matrices S± and T± are defined in Eqs. 
and H2.11|l and we introduce the matrix function 

T±{Na,m) = ^ ^-('+1)5-^1 (; + i)i?,*± (/)£;', 



^±iz) = T±(-c»,oo). 



(6.5) 



Note that they are the analytical solutions 'i'± that enter 
naturally the matrices T±. 

We have 



Now we derive a variation of ^' 



+ • 



from Eq. (ESJ that 5*+ = SJ+E'^S+E-"" + 
J+E"5S+E^''\ The first term in r.h.s. is transformed 



to ie'^+{n)E''T+{n,oo)E-"St by means of Eq. (j^ . 
while the second term, due to Eq. (|6.4|l and a trick 
(55+ = (55* Mil, Mil = diag(l,0), is written as 
-ie^+{n)E''T+{z)MiiE-''6t. Therefore, the variation 
of ^'+(n) takes the form 



St 



-i£^+{n,z)E'^Il+{n,z)E-"-, (6.6) 



where H_|_ is the evolution functional [3q| defined here by 



, , / T+ii(-oo, n - 1) -T+i2(n, oo) 
z) - ^ x+2i(-oo, n - 1) -T+22(n, oo) 

. (^-^^ 

Therefore, in the case of perturbations the evolutionary 
equation for 4"+ gains the additional term responsible for 
the perturbation: 



^r+i = V'i+ - - ie'i'+E'm+E-'^ 
Similarly, 



St 



ie£:"H_£;-"«'" 



(6.8) 



(6.9) 



with 



H_(n, 2) = 



and 



T_ii(-cx),rt - 1) T_i2(-oo, n - 1) 
-T_2l(?^, 00) -T_22(", 00) 



(6.10) 



^-ij = V + f^fl^ + ie£;"n_£:""*l\ (6.11) 

Remarkably, the functions T± are interrelated by the 
matrix G entering the RH problem: 



GT+G- 



(6.12) 



The evolution functionals Il± play the key role in the 
analysis of a perturbation because they contain all needed 
information about it 's^. It is seen from the defini- 
tions H6.5|l . (|6.7|l and l|6.10|l that the matrices H± are 
meromorphic (and E"Il±E~" are bounded) in C± hav- 
ing simple poles at zeros of det^+(z) and det^"!^, re- 
spectively. Further, the evolution equation for G follows 
easily from Eqs. (|3.1|l . (|6.8I) and (|6.11() and takes the form 



Gt = [n,G]~ie{GU+-Il^G), (6.13) 
or, for G — exp(— f7i)G'exp(f7i), 

Gt = -ie(Ge""*H+e"* - e-"*n_e"*G). (6.14) 

B. Discrete data 



In the point zi 



^+{n,zi)\n) = 



(6.15) 
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and near this point 



VII. ADIABATIC APPROXIMATION 



where II^''*'' is the regular part of 11+ in the point zi and 
ReSz=2i stands for the residue at z = zi. It is shown 
in the Appendix A that the evolution equation for the 
eigenvector takes the form 

= f](zi)|n)+zeS"(zi)nV^^^)(zi)^;-"(zi)|n}. (6.17) 

Remember that the n-dependence of \n) is given by 
Eq. H4.10|l . \n) = E'"-{zi)\p), with the n-independent vec- 
tor \p). Therefore, the perturbation- induced evolution 
of the vector \p) = exp[—n{zi)t]\p) is governed by the 
equation 



z — Zl 



-Res,=,,n+(z), (6.16) 



(6.18) 



In the absence of perturbation, the vector \p) in 
Eq. I|6.18|l coincides with that in Eq. 14.11|l . 

Evolution of zero zi is derived by taking the total time 
derivative of det^'+(zi) — 0. We obtain 



zit 



jd/dt) det^+(z) 
{d/dz)det^+{z) 



Because (5/9z)det*+ = -ietrn+ det ^'+, dot = 
v+{n)a+{z) [see Eq. (|TTUI) ] and a+{z) = (z^ - z1){z'^ - 
Zi)~^ , the latter formula following from a+(z) = a+(— z), 
lima(z) 1 as z ^ c» and a+(±zi) = 0, we obtain a 
simple equation 



zit = zeReSz=zitrn+(7i, z). 



(6.19) 



It is important that l.h.s. of Eqs. and (|^?T!?|) 

do not depend on n. Therefore, we can consider these 
equations for n +00 where 



n+(z) = 



T+ii(z) 
T+2i(z) 



Within the adiabatic approximation, we ignore radia- 
tion effects and assume that the soliton adjusts its hy- 
perbolic secant shape to perturbation at the cost of slow 
evolution of the parameters. Evolution equations for the 
soliton parameters in the adiabatic approximation have 
the form 

/i( = e sinh /i (7.1) 
Im(i?„) cosh /i(rt — x) 
cosh /i(n + 1 — x) cosh /i(n — 1 — a;) ' 

kt — —tsmh^ (7.2) 
Re(i?„) sinh/x(n — a;) 



E 

1 — — cx 

00 

E 

n— — oc> 



cosh /i(n + 1 — x) cosh /i(n — 1 ~ x) 



2 . . e . 

Xt = — sinh /i sin fc H sinh /i 



(7.3) 



E 



(ri — a;)Im(i?„) cosh /i(n — x) 
cosh ^{n + 1 — 2;) cosh ^{n — 1 ~ x) 



at = 2(cosh ^ cos fc H — sinh ji sin fc — 1) (7.4) 

C30 

-|- e I [(n — x) sinh /i sinh/i(n — x) 

Tl — — C30 

— cosh^cosh^(n — a;)]Re(i?„) 

A; ~| 
+ — (n — x) sinh/iCOsh/i(n — .T)Im(i?„) > 

X sech^(n -I- 1 — a:;)sech^(n — \ — x). 



Here _R„ = r„ exp[— ifc(n — x) — ia] and r„ is constructed 
by means of the AL soliton solution (|5.9(l . Eqs. (|7.1II- 
(|7.3|) have been obtained for the first time in Ref. 2;^. 
The derivation of Eqs. (|7.1|I - H7.4|I within the RH problem 
approach is given in the Appendix B. 



As a result, the evolution equations for the discrete RH 
data are finally written as 



zit 



\p)t 



ieRes,=,,T+n(^), (6.20) 
Xzi)t ( "^+1? (^1) 1 „^(^i)t\ 



zee 



T^'^f(zi) 



'b).(6.21) 



It should be noted that Eqs. (EHI), (|On|) and ^^1^ 
are exact. However, they cannot be directly applied be- 
cause the matrices H± and T± depend on unknown so- 
lutions ± of the spectral problem with the perturbed 
potential. In the following sections we will describe for 
sufficiently small e the iterative RH problem-based pro- 
cedure to consecutively account for two main approxima- 
tions: the leading-order adiabatic approximation and the 
next-order (the first-order) one. 



VIII. RADIATION EFFECTS 

The continuous part of the RH data describes a distor- 
tion of the soliton shape and emission of small-amplitude 
dispersive waves by soliton. To account for the continu- 
ous data, we should abandon the condition G = 1 and 
admit a z-dependence of the regular RH problem solu- 
tions i]-)±. In other words, we pose 

G = 1 + eg{z), V'+(", z) = ^+(«)(1 + £'/'(", 2)), 

(8.1) 

where -0" stands for the solution H5.1|l of the regular RH 
problem (|4.3|l in the adiabatic approximation, whereas 
the off-diagonal matrices g{z) and (j){z) describe first- 
order corrections. Therefore, the reconstruction for- 
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mula H4.9(l takes now the form 



lim 



r[\^ (n)r22 (n, 0) - ecj^iy H^22 (n, 0) 



[zV'+(l + e0)r] 
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(8.2) 



(1)/ 



The first term in the r.h.s. of Eq. (|8.2(l represents the 
famihar sohton solution in the adiabatic approximation 
and the second one is responsible for radiation (soliton 
shape distortion). For the derivation of Eq. (|8.2|) we cm- 
ploy the fact that the off-diagonal matrix satisfies the 
asymptotic condition z^^^^^^ -I- 0{z~^) with 



Evaluation of <p^^2 hence of radiation correc- 

tions to soliton solution reduces to solving the regu- 
lar RH problem ijHSll with G as in Eq. (|8.1|l . In- 





^(1) 



deed, we have 



1 -I- eri;"g(z)£;-"r-i and 



the jump of the piece-wise holomorphic function ipi^) = 
{ip+{z), z £ C^;tjj_{z), z € C-} across the contour \z\ — l 
is written as 



(8.3) 



Here we omit terms with higher order of e and invoke the 
equality ^p'^ = [see Eq. (|5.1|l ] valid in the adiabatic 
approximation. The Plemelj formula gives for z G C+: 



1 + 



_1 z' — z 



(r£;"5£;-"r-i)(z') 



Inserting here ip^ from Eq. (|8.1|) and performing the 
asymptotic expansion at z — > cxd, we obtain the expan- 
sion coefficient 



i)(n) = — — I dz(r£;"g£;-"r-i)(z) 



(8.4) 



determining the radiation correction H8.2(l . Therefore, 
our next step is concerned with finding the matrix g. 

To this end, we turn to the evolution equation (|6.14|) 
for the matrix G which is evidently related to g: 



G = 1 + eg, 



(8.5) 



Substituting this equation into Eq. (|6.14() . we obtain in 
the first order of e 



Moreover, it follows from Eqs. (|6.5|l and (|8.1() that T_ 
T+ in the first order of e. As a result, 



-nt 







-^ + 12 





and the equation for gi2 takes the form 

5i2t = i exp \~i{z - z^'^ft] T+i 



■7) 



It is important to stress that because T+12 corresponds 
to the first order correction, we can replace in the defini- 
tion 1)6. 5|l of T-|_ unknown solution ^+ of the regular RH 
problem 1)8. 3|l by the known one if)]^. Integrating then 
Eq. 1)8. 7|l . we can find the matrix g (|8.5|) . 

The further stage is to consider the integrand in 
Eq. ||S3Il. It can be shown from {TE'^gE^''T-^)i2 = 
z-2'Ti2(r-i)i2.g2i + z2"rii(r-i)225i2 and explicit ex- 
pressions 1)5. 7|l for r that the term with 521 is multiplied 
by sech^ /i(n — a; — 1) and hence vanishes at n ^ ±00. As 
a result, we are left with 



/12 = (rs"5S-'T-i)i2 = 



^2n „ 

512, 



1 ^2n „ 
t2 512, 



(8.8) 

Let us summarize the main steps in calculating the 
radiation correction for a given perturbation r„. First, 
we should explicitly find the function T_|_i2(z) from the 
definition H().5|) with — ip^T, V'+ and F being given 
in Eqs. H5.1|l and H5.7|l . respectively. Then we integrate 
Eq. H8.7|) and obtain the matrix g given in Eqs. (|8.5|) and 
(|8.1|) . For the known function 512 (2^) we obtain the inte- 
grand (|8.8|l . Finally, after calculating the integral H8.4I) 
we arrive at the needed result. 

In the next section we illustrate the proposed formal- 
ism on an example of calculating the radiation corrections 
to the AL soliton in the case of some model perturbations. 



IX. EXAMPLES 

Here we apply our formalism to describe the perturbed 
AL soliton dynamics for the typical representatives of dis- 
sipative and conservative perturbations - linear damping 
r„ = —iun and quintic perturbation r„ = |u„|^u„. The 
interplay between the dissipative and conservative per- 
turbations for the AL model is considered in the adiabatic 
approximation by Abdullaev et al. j29j | and numerically 
by Soto-Crespo et al. [iH. 



-Sit 



(H+-n_)e 



nt 



(8.6) 



Because g does not depend on n, we can put n ^ 00 in 
Eq. 1)8.6)1 which gives 



n+(n^oo)-H_(n^oo)= ( ^"^^ 

' JL +21 tl 



A. Linear damping 

In this case Rei?„ = 0, Imi?„ = — sinh /i sech/i(?T, — x) 
and we have in the adiabatic approximation 



fc = const, sinh/i = sinh(/io)e 



-2et 



Mo = = 0), 



n ' 1 1 1 ' 1 ' 1 

5 10 15 20 t 

FIG. 2: Evolution of the soliton mass (m — 2fi) and group 
velocity Vgr in the case of linear damping for k — 7r/4, e = 0.03 
and fi{t = 0) = 1. 



.5 5 10 15 20 n 

FIG. 3: Evolution of the perturbed AL soliton u„ 18.211 with 
account for the first-order correction. The soliton parameters 
are the same as in Fig. 2. 



Xt 



Ugr 



2-Ke tanh 
/i^ sinli(7r^//x) 



sin27rx, Sa/St — kxt- 



In the process of obtaining the equation for xt we use the 
Poisson summation formula |40j 



E = - / dyfiy) 



cos ■ 



2TTsy 



1-1 ^^-^^ 
and, following Ref. |2J|, we restrict ourselves to the lin- 
ear harmonic term (s = 1) only. Higher harmonics con- 
tain the factor exp(— 7r^s//x) which for /i « 1 is evidently 
small. Hence, mass of the soliton decreases exponentially, 
its group velocity acquires a constant value (= 2sinfc) 
after some transient period (Fig. 2), while its phase is 
governed by the evolution of the soliton position x{t). 

Now we embark on a calculation of radiation effects. 
Following the prescriptions of Section IVIIII we find at 
first the matrix function T+ written in accordance with 
Eq. Ht).5|l as 

oo 

T+(z)= E £^~^"+^H« + l)^"^(r^+l)7^„^(r^)£;". 

n— — OO 

(9.2) 



Here 



and 



Tin = (V'°)-i(n + l)i?V'[(n) 

r„r22(n,0)- 
<r22(n-|-l,0) 

r22(n, 0) = 



cosh /i(n — x) 

Substituting F and F^^ (|5.()|) into Eq. H9.2|l . we arrive at 

^ — fi+ia — ikx 



T 



-12 



cosh/x — cos(fc — 29) 
(l — cosh//cos(A: — 20)) S*! + 2 sinh sin(/c — 26')S'2 



where 



OO 

E 

n— — OO 



2n ( cosh ii{n—x) \ 
1 sinh fi(n — x) J 



cosh /i(n — X — 1) cosh fi(n — a; + 1) 



Calculating these sums by means of the Poisson formula 
(|9.1|l . we obtain a simple expression 



T 



+ 12 



exp(— /J, - 



2i0x) 



fiz cosh fi cosh(7r(fc — 26) /2^) 



.3) 



Here we pose z — exp(i0) bearing in mind subsequent 
integration along the contour \z\ = 1. What is more, 
because the radiation correction (|8.2I) is multiplied by e, 
we restrict ourselves to the leading term in each sum. 
Integrating then Eq. H8.7(l for gi2 with T+12 of the form 
Eq. H9.3|) and transforming the result to 1712 in accordance 
with Eq. (|8.5() . we get 



512 



exp(ia — 2i9x) 



1 _ e-*A(e)t 



cosh cosh(7r(fc - 26) /2^) A{9) 



Here A{6) — {k — 29)vgr + 2(cosh/icosfc — cos26') and, 
within the first-order approximation, we can take as Vph 
and Vgr their initial values. Therefore, we arrive at the 
integrand I12 (18.811 which determines the integral (|8.4|l . 
This integral can be calculated by residues. The domi- 
nant contribution is provided by the third-order residue 
in the point z = zi (note that both cosh(7r(fc — 26) /2^) 
and A(6) have simple zero in this point). The result- 
ing expression is rather lengthy and does not reproduced 
here. Instead we plot in Fig. 3 the evolution of the 
perturbed AL soliton w„ (|8.2|l with account for this ex- 
pression. The shape of the soliton changes periodically 
due to the discreteness of the system, with simultaneous 
decreasing of mass in virtue of damping. More detailed 
conclusion about the first-order contribution can be in- 
ferred from Fig. 4 where a difference |u„| — \us\ is pic- 
tured. Here Ug is the AL soliton in the adiabatic approx- 
imation. The shape distortion is mainly localized within 
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|u„Hu.l 




FIG. 4: The difference \u„ \ — \us\ for different time intervals 
indicative of the shape distortion effect. Here Us is the sohton 
in the adiabatic approximation and the sohton parameters are 
the same as in Fig. 2. There are no long-lived nonvanishing 
dispersive waves. 



Radiative corrections for the quintic perturbation are 
much the same as for damping. Indeed, the function 
T_(-i2 takes the form 



smhfi 



cosh(7r(fc - 20) /2n) 



where 



'3(0) =i (cosh A* 



sinhi 



(1 - I sinh^ ^) (cosh ^ + sinh A*) + f cosh / 



- (^2 + -sinh> 



the soliton 'envelope' and in general is asymmetric. Such 
a behavior agrees with the asymptotic representation of 
the first-order correction for n — > +00: 



'«n(rad)=-e0l2^r22(ri,O) 



-—tanhne-^'+"'+'''{n-x) 
2H 



sin k 



1 



sinh 2fi 



H- 



sinh/i y 2/i 

(9.4) 

Here H — sinh /i cos k — i cosh /i sin k + sinh /i sin k. 
It follows from this expression that there are no nonvan- 
ishing linear waves at n ^ 00. The similar result takes 
place for n ^ — 00. 



B. Quintic Perturbation 



k-26 



■ sinh / 



-(1 + cosh^ /i) — cosh /.t cos(fc — 29), 



Pie) = [sinh^(cosh/x - cos(fc - 2e))]~^ . 
Therefore, the integrand I12 is written as 

J ^ Sinh^ ^^-^1 2(»-x)-l 

2/x cosh(7r(fc - 29) /2^i) ~ zf 

P3(^) -2e*('=-29) 



A{9) 



with the same A(6'), as before. The result of calculation 
of the integral (|8.4(l with the above integrand has the 
same structure as in Eq. (|9.4|) . What is more, the func- 
tion P3{9) — 2e*'-*'~^^'' being zero for z = zi, does not 
contribute to the n^-order of the radiative correction. 



It follows from the results of the Section IVIII that for 



Re(i?„) — sinh^/i sech^/i(n — x), Im(i?„ 
we have in the adiabatic approximation 



/i = const, 

2e7r sinh* fi 
-ifJ- sinh — 



sinh^ . , 



M* 



27^ 
/i2 



sinh /X 



kt = (9.5) 

— sin27ra;. 
3 



We do not write here the evolution equation for the phase 
a because of its lengthy form, though it is found quite 
immediately. Hence, soliton mass is preserved in the 
presence of the quintic perturbation, while the param- 
eter k{t) (and hence the group velocity) oscillates with 
a very small amplitude near the initial value. Linear 
stability analysis demonstrates that the fixed points of 
the evolutions (|9.5() kg = ttto, m = 0, ±1, ±2, . . ., and 
Xg = 1/2, I = 0, ±1, ±2, . . . are stable for even (odd) m 
and odd (even) I. 



X. CONCLUSION 

We have proposed a formalism suitable for analytical 
investigation of dynamics of the AL soliton subjected to 
a perturbation. This formalism provides a possibility of 
calculating both evolution of the soliton parameters and 
perturbation-induced radiation effects. Remarkably, it is 
the RH problem-based approach that has been proved to 
be efficient for treating continuous nonlinear equations, 
both integrable and nearly integrable, that turns out to 
be the natural basis to study discrete nonlinear systems. 
We have demonstrated within this approach how to con- 
sistently advance from an integrable to perturbed system, 
in so doing the only ingredient that should be added to 
the formalism to account for a perturbation is the evolu- 
tion functional H+ (or n_) introduced by Shchesnovich 
|36| . A natural step to further extend the applicability 
of analytical methods in the the theory of discrete non- 
linear systems is to consider vector AL-type solitons 42| . 
Work in this direction is now in progress. 
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APPENDIX A: PERTURBED EVOLUTION OF 
EIGENVECTOR 

In this Appendix we derive Eq. (|6.17() . Taking the 
total time derivative of Eq. H6.15|l gives with account of 
Eqs. (ESI) and 



(z-zi)-iRes,^,,(i;"n+£;-") 

d 



zt — \n) + ^+{n,zi)\n)t^O. (Al) 

OZ J z=zi 

Let us introduce a holomorphic function 11 = — ie(z — 
Zi)E^Il^E~" which evidently gives 

fl(zi)|n} = -ieRes^=^, [E"U+{z)E^"-] \n). (A2) 

On the other hand, representing 11+ (z) from Eq. H6.8|l as 

-ieE''U+{z)E'" = I'+^'J'+t - + n, 

we obtain 

n(zi)|n) = [iz ~ zi)*;i*+4^ \n) = -zu\n). (A3) 

Comparing Eqs. ljA2p and (|A3p . we arrive at 

ieRes,=., [i;"n+(z)S-"] |n) = zit\n). (A4) 

Applying now ^"+(21) to both sides of Eq. ljA4|l . we ob- 
tain the important identity 

^+{n, zi)Res,=,, {E"n+E-'') = 0. (A5) 



Eqs. (jA4|) and (|A5|) permit to considerably simplify 
Eq. (|Aip . Indeed, the last term in square brackets in 
Eq. ijHJ is rearranged by means of Eq. ljA5|l as 



-zit 



z — Zl 



Res,=.,(i;"n+(z)i;-") 



\n) 



Res,=,,(i;"n+(z)i;-")|n) 



\n) 



and cancels the same term in Eq. (|A1|I . As a result, the 
evolution equation for the vector \n) takes the form 

\n)t = n{zi)\n)+ieE'^{zi)U^l''^\zi)E~^{zi)\n). 



Here we obtain within the RH problem approach 
Eqs. H7.1|l - H7.4|l which govern the adiabatic dynamics of 
the AL soliton. 

First of all we turn to Eq. (|6.20|l . In accordance with 
Eqs. I|6.5(l . H4.1|l . (|5.1|l and (|5.2|l we write 

Res2=2iT+ii(z) = Res^^^i 

- J2 ^'\n+l,z)^^\n+l)R„^+in)T{n,z) 



11 



sinh^ 
2zi 





r22(71+l,0)< 



n— — 00 







r(n, Zl) 



With account of explicit expressions H5.6|l and l|5.7|l for 
F- and F we obtain 



zit ~ ^—zi sinh/i 



E 



n— — oo 



cosh + 1 — .t) cosh — 1 — x) 



where i?„ = r„ exp[— ifc(n — x) ~ ia] . Then from the defi- 
nition Zl = exp[(l/2)(/i + ifc)] we easily derive Eqs. 1)7. 1|) 
and ifT^ . 

In order to obtain Eqs. (|7.3(l and (|7.4|1 . we should at 
first calculate T("s)(zi): 

T("s)(zi) =[T(z) - (z - zi)-iRes,^,,T(z)]^^ (Bl) 

sinh /i 



smh/i ■r-^ 
lim — r > 



^ F+(n-)-l) 7e„F(n,zi)£;" 
4zi / 



n— — 00 



z^zi 2(z — Zl) 

X F(n,z)£;"(z) -S-("+i)(zi)F_(n + l)7e„F(n,Zi) 



Here 7?.„ = ip^Hn + l)^n'0+("')- The second term in the 
r.h.s. of Eq. IjBip gives the ratio 0/0 in the limit z zi. 
Using the I'Hospital rule, we ultimately arrive at 



Y(rcg)^ 



OC 

E E-i-+'Hzi)l{t 

1.— — nr> ^ 



sinh// , , sinh// , 
—i^a^F^in + 1 + — n + 1 7^„^ n, zi 
2zi 4zi 



sinh II 
+ ^ 

sinh u 

— n 

2zi 

Therefore, 



F-(n) 



F+{n) 



^(zi-zi)2 (zi-fzi)2 
[CT3,i^-(n + l)7e«r(n,zi)]ls"(zi). 
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(3i?„e^("-^) - iJ^e-'^^"-^)) sinh/i + 2(i?„ cosh/x + Ke-^')e 



-4i?„sech^(n + 1 — x) 



sech/i(ri + 1 — a::)sech/i(n — 1 — x), 



(B2) 



4^^^^cosh/i(7T, + 1 — a;) cosh /i(n — 1 — a;) 



{ cosh M ^ 2 ^i'^^ M j 



(B3) 



+ 2e-^("-i-^) cosh^(n - x) + ie-2M(n-x) sinh^)i?; - 2risinh^ (i?„ + e-^'''^"-^'^ R*„ 



Then we obtain from Eq. H6.21|l the following evolution 
equations for the components of the vector \p): 

Pit = ieT^\°'^-'(zi)pi, 



P2t = ieT^\°^^(zi)exp 



i J {zf + Zj"2 - 2)dt 



Pi 



Because {pi/p2) = exp(a + itp), we have from Eq. (jB2|) 
d 



(a + iip) 



" 7 X] (3i?„e^("-^) -i?;e-''("-^)) sinh/x 

+ 2n (i?„e''("-^) - i?;e-''("~=")) sinh/x 

- 2i?* cosh fi{n - a;) 

X sech/x(n + 1 — x)sechfi{n — 1 — x) 
which results in 



"+2 



at = — esinh/i 

n— — c 

Im(_R„) cosh — a;) 
cosh/i(n + 1 — a;) cosh/i(ri — 1 — x) 



(B4) 



Inserting here Eqs. (|7.1|l . H7.2|l . (|B3I) and HB4p . we finally 
obtain Eqs. (|73|) and (|T^ . 



(ft = ^ sech/i sech/i(n — x) — cosh /i cosh /i(n — a;) 



n— — oo 



+ — sinh/x sinh /i(n — a;)] Re(i?„) 
X sech/i(n + 1 — x) sech(7i — 1 — a;). 
It follows from Eqs. I|5.8|l and (|5.1U|) that 



. ^^inh,sinfc-i 



^+2 ) A^t + a* 



at = 2 I cosh /i cos fc H — sinh /i sin fc — 1 

3\ fc fc 



x + -]kt 



^ + 7; -A^t + ft- 

2J fi 
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